Two super Tonks- Girardeau states of a trapped ID spinor Fermi gas 
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A harmonically trapped ultracold ID spinor Fermi gas with a strongly attractive ID even-wave 
interaction induced by a 3D Feshbach resonance is studied. It is shown that it has two different 
super Tonks- Girardeau (sTG) energy eigenstates which are metastable against collapse in spite of 
the strong attraction, due to their close connection with ID hard sphere Bose gases which are highly 
excited gas-like states. One of these sTG states is a hybrid between an sTG gas with strong (^4-) 
attractions and an ideal Fermi gas with no ("ft) or (44) interactions, the sTG component being an 
exact analog of the recently observed sTG state of a ID ultracold Bose gas. It should be possible 
to create it experimentally by a sudden switch of the (t4-) interaction from strongly repulsive to 
strongly attractive, as in the recent Innsbruck experiment on the bosonic sTG gas. The other is a 
trapped analog of a recently predicted sTG state which is an ultracold gas of strongly bound (t4-) 
fermion dimers which behave as bosons with a strongly attractive boson-boson interaction leading 
to sTG behavior. It is proved that the probability of a transition from the ground state for strongly 
repulsive interaction to this dimer state under a sudden switch from strongly repulsive to strongly 
attractive interaction is -C 1, contrary to a previous suggestion. 
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If an ultracold atomic vapor is confined in a de Broglie 
wave guide with transverse trapping so tight and tem- 
perature so low that the transverse vibrational excitation 
quantum is larger than available longitudinal zero point 
and thermal energies, the effective dynamics becomes 
one-dimensional (ID) [l], 3D Feshbach resonances 
Q allow tuning to the neighborhood of ID confinement- 
induced resonances 0, H| where the ID interaction is 
very strong, leading to strong short-range correlations, 
breakdown of effective-field theories, and emergence of 
highly-correlated TV-body ground states. In the bosonic 
case with zero-range repulsion gs5(xj — xg) with coupling 
constant qb — > +oo, the Tonks-Girardcau (TG) gas, the 
exact iV-body ground state was determined in 1960 by a 
Fermi-Bose (FB) mapping to an ideal Fermi gas lead- 
ing to "fermionization" of many properties of this Bose 
system, as recently confirmed experimentally [HQ- It 
is now known [8-10] that the FB mapping is of much 
greater generality; when supplemented by an inversion 
and sign change of the coupling constant, it provides a 
mapping between the iV-body energy eigenstates of a ID 
Bose gas with delta-function interactions gsS(xj — xg) of 
any strength [Lieb-Liniger (LL) gas [Tl|] and those of a 
spin- aligned Fermi gas. 

A generalization of the FB mapping to a ID spinor 
Fermi gas is relevant to a recent prediction of a super 
Tonks- Girardeau (sTG) state in such a system [l2j], and 
will be employed herein. The model consists of an ul- 
tracold gas of fermionic atoms in an effectively ID trap, 
in two different hyperfine states which are conveniently 
labelled as t (spin up) and 4- (spin down). The previ- 
ous analysis [12| assumed a ring trap (periodic boundary 
conditions), but all experiments on ultracold gases with 



tight transverse trapping (ID regime) use a linear geom- 
etry with a longitudinal harmonic trap potential. The 
Hamiltonian generalizes that of jl2| by adding a spin- 
independent harmonic trapping potential: 
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(1) 

Since spin-up atoms are distinguishable from spin-down 
atoms, both 3D s-wave and p-wave scattering are allowed. 
The case relevant to generation of a sTG state is that of 
strong s-wave scattering due to a 3D s-wave Feshbach 
resonance, which leads in ID to a confinement-induced 
even-wave resonance [H, HJ and a LL delta function in- 
teraction [Til] - It acts between atoms j and £ only if 
the wave function is symmetric under spatial exchange 
(xj <-> xi), in which case fermionic antisymmetry un- 
der combined space-spin exchange (xj, Oj) -H> (xi, at) de- 
mands that the spins be antiparallel, Uj =f and <jg =\, or 
vice versa. On the other hand, the delta function interac- 
tion is cancelled in the case that the spin orientations are 
tt or 44-i in which case the required spatial antisymmetry 
requires that the wave function vanish when Xj — xg. 

If gp is large and negative, two quite different sTG 
states can occur. One is a hybrid between an sTG gas 
with strong (fi) attractions and an ideal Fermi gas with 
no (tt) or (11) interactions, the sTG component being 
an exact analog of the sTG state in an ultracold ID Bose 
gas predicted in (TH-HH and recently created by the Inns- 
bruck group [lj|. It should be possible to create it ex- 
perimentally by a sudden switch of the (ft) interaction 
from strongly repulsive to strongly attractive, as in [llij . 
The other sTG state is a trapped analog of a recently pre- 
dicted |12| sTG state which is an ultracold gas of strongly 
bound (t|) fermion dimers which behave as bosons with 
a strongly attractive boson-boson interaction leading to 
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sTG behavior. Both of these sTG states will be studied 
herein. 



Consider first the case N = 2 
parallels that for the spinless Bose gas 
interaction is gp8{x\ — X2), the harmonic trap poten- 
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Huj/2 assuming that the cm. 
unexcited. The form of the relative wave function </)(x), 
with x = Xi — X2, now depends on whether the spins o\ 
and o"2 are parallel or antiparallel. In the parallel case the 
relative wave function <j>(x) is antisymmetric so that the 
LL interaction is cancelled, and </> reduces to a trapped 
harmonic oscillator eigenstate odd in x, the product of 
the ground state Gaussian and a Hermitc polynomial of 
odd order. In the antiparallel case <fi(x) is symmetric and 
is identical with the solution for spinless bosons pjHUj]. 
The solutions are analytic continuations of the Hermitc- 
Gaussians to nonintegral quantum number u, and are 
parabolic cylinder function D u (q) where x = qx osc and 
the allowed values of v are solutions of a transcenden- 
tal equation F(h — \v)/Y(—^v) — — A where A is a di- 
mensionless coupling constant A = gp /2 3 / 2 hu)x osc [l8j]. 
This gives the solution for q > 0, while for q < it is 
D u (\q\), since 4>{x) is even. For Xp — » +00 the solu- 
tion reduces to the first excited harmonic oscillator state 
<j)(x) = \q\e~ q I" 2 [2(j, the usual TG gas ground state with 
a cusp at q = due to the point hard core interaction. 
For A — > — 00 this wave function is still an exact energy 
eigenstate, but it is highly excited, the ground state be- 
ing a collapsed state which is an analog, for the trapped 
system, of McGuire's cluster state (22[. It is an even so- 
lution also expressible in terms of a D v , but one whose 
energy approaches —00 as gp — > —00 (<zitj — > 0+); see 
Fig. 5 of [lg]. For aiD — > 0+ it is well approximated by 
« exp(— \xi - x 2 \/an)) exp[-(a;^ + x\) /2x 2 os i\ where 
> is the ID scattering length. aiD vanishes as 
-> 00 and increases with decreasing |A| [l|, 0, [l9j . 
at the same time the node moves from the origin 
to a position x no d which is close to am when |A| 3> 1 
but slightly smaller, the explicit expression being given in 
! 17t| - A s m the case of spinless bosons, the wave function 
for I a; I > lEnode is identical, apart from normalization, 
with that of a ID hard sphere Bose gas [5j with hard 
sphere diameter rr no de [171 • 
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by wave functions of the form 

Ipp = M(xi,ai\"- ;X A r,(TAr)V'ideal 
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where the spin-dependent Fermi-Fermi mapping function 
M, which maps the strongly interacting spinor Fermi gas 
to a spinless ideal Fermi gas, is 

M(x 1 ,a 1 ;---;x N ,a N )= JJ a(Xj,aj;Xi,ae) 

i<j<e<N 

a(xj,aj;x£,ai) = (S ajt 5 ae i - 8 a]i 8 af/[ )sgn(xj - xi) 

+8<r i t8<r t t + So-j^o-tl , (3) 

the signum function sgn(a;) is +1 ( — 1) if x > (x < 0), 
and Videai is an energy eigenstate of the trapped ID 
ideal gas of "spinless" fermions, a Slater determinant of 
N different harmonic oscillator orbitals. Starting from 
this ideal Fermi gas, this mapping generates the required 
TG-gas cusps in the singlet (spin-antisymmetric, space- 
symmetric) channel of the strongly interacting spinor 
Fermi gas, but no interaction cusps in the triplet (spin- 
symmetric, space-antisymmetric) channels. Here we are 
particularly interested in the ground state i/jfq, which is 
mapped from the ground state V'ideaio of the trapped ideal 
gas, a Slater determinant of the lowest N single-particle 
eigenfunctions 4> n of the harmonic oscillator (HO): 



idcalO 



(xi, ■ ■ ■ , X N ) = 



1 



(N-1,N) 

det 4> n (xj) 



'N\ («.i)=(o,i) 



The HO orbitals are 
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with H n (Q) the Hermite polynomials and Q = x/x osc , 
where x osc = y^h/mui is the harmonic oscillator width. 
The ground state is a van der Monde determinant [2(| 

(N-1,N) (N-1,N) 

det H n { X j) = 2 N ^ N -^' 2 det ( Xj ) n 

(n,j)=(0,l) (n,j) = (0,iy 

l<3<t<N 

yielding an exact analytical expression of Bijl-Jastrow 
pair product form for the A^-fermion ground state: 



ipFo(xi,cxi', ■ ■ ■ ;x N ,a N ) = C 
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x Yl a(xj,aj;xe,ai)(xj - xi) 
i<j<e<N 



Generalizing to N > 2, denote the iV-atom wave func- 
tions by tPf(xi, o"i ; • • • ; xn, cn)- Consider first the case 
A — > +00. In this limit the wave functions are required 
to vanish if Xj — xi provided that <jj =t , ere —i or 
vice versa, and if cr j = at they vanish automatically at with normalization constant 
Xj = Xf. by antisymmetry. Furthermore, if the spatial 
coordinates {x\, ■ ■ ■ ,xn) arc all different, then the inter- 
action vanishes independently of (a\,--- ,ctn) and the 
wave functions are kinetic energy eigenstates which are 
symmetric under Xj Xn exchange if crj 7^ ae but anti- 
symmetric if crj = <T£. These requirements are all satisfied 
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This TG-ideal Fermi gas hybrid is the exact ground state 
in the TG limit A = +00. 
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Suppose that now A is changed instantaneously to — oo. 
Then the state ([7]) is still an exact energy eigenstate, since 
H commutes with the mapping function M except at 
unlike-spin collision points xj = Xg, where the wave func- 
tion vanishes, and the state ([7]) maps to the ideal Fermi 
gas d, H[ as in the bosonic case d, H, [l7|. In fact, this is 
an exact energy eigenstate in the limit |A| — ¥ oo even in 
the dissipative case where A is complex, since the wave 
function vanishes at contact in that limit [21j , validating 
the FB mapping. This follows from the LL contact condi- 
tion pi[ 2[d^/dxjk] Xjk =o+ = (m.9F/ft 2 )V'(0) with a cusp 
(derivative sign change) at the origin. Assuming tp nor- 
malized so that [dtp/dxj k ]x Jk =a+ = - [dip/dx jk ]x jk =o- = 
1, this leads to the Taylor expansion ip = (2H 2 /rngp) + 
\xik \ + • • • about X j — X y so ip vanishes at contact when 
qf — > oo. For A -C — 1 such a state is highly excited. The 
much lower ground state is a hybrid of a totally collapsed 
McGuire cluster state [22j for antiparallel spins and an 
ideal Fermi gas ground state for parallel spins, and Tay- 
lor expansion fails in the limit A — > — oo (<2id 0), 
where the derivative of e~> Xi ~ Xk '' aiD diverges. Although 
the exact ground state is not known for N > 2 and 
— oo < A <C — 1, one expects that it will have the form 
of Eq. (3) of [17| for pairs with antiparallel spins, while 
for pairs with parallel spins it will be an ideal Fermi gas 
ground state. One therefore expects that it will be well 
approximated by 
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(9) 



In fact, the LL contact conditions are satisfied exactly at 
each collision point Xj = X£. 

Suppose that A is switched from large positive to large 
negative values rapidly enough that the sudden approx- 
imation is valid. Then the initial wave function will be 
nearly equal to the state of Eqs. and ([7]), and the 
wave function after the switch will be a superposition of 
all eigenstates of the system with the given negative A 
value, Mt) = T, a {^xMF^x a e- lE ^ t ' h . The domi- 
nant term in ip\ will be an sTG state which reduces to 
(UJ as A — s- — oo. The obvious generalization of the N=2 
sTG wave function to arbitrary N differs from Eqs. @ 
and (J7J through replacement of the ideal Fermi gas factor 
(xj — xi) by the N = 2 solution D u for antiparallel spins, 
while retaining the ideal Fermi gas form for parallel spins: 
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where 

+(^ it <W ~ S^ tt )D v (\q je \)e^ (11) 

and qji = (xj — xi)/x osc . It satisfies the contact con- 
ditions exactly. One expects the existence of a highly 
excited gas-like sTG state with nodes only at a nearest 
neighbor separation \xji\ = 1^ when the spins are an- 
tiparallel, as in the Bose case [l7|, where x no dc increases 
with decreasing |A|, is very close to am for |A| 3> 1, and 
goes to zero along with am in the TG limit |A| — > oo. 
For all N > 2 the approximate wave functions (fT0|) have 
exactly these properties, vanishing only at \xji\ = x no d c 
when the spins are antiparallel, and becoming exact both 
at the collision points Xjt — and when all \xji\ — > oo. 
Hence we expect the unknown exact sTG excited state 
for finite negative A to be well approximated by the state 
ofEq. fTO]). 

To establish a connection with a hard sphere ground 
state, consider first the case N = 2. If both atoms have 
spin up or both have spin down, then the wave functions 
are necessarily spatially antisymmetric, vanishing when 
x\ = x%, so the delta interaction in Eq. (JTJ) is ineffective 
and the energy eigenstates are those of the trapped ideal 
Fermi gas. On the other hand, if they are antiparallel 
in a spin-antisymmetric singlet state arising from an s- 
wave Feshbach resonance, then the ground state wave 
function is spatially symmetric and of the previously- 
discussed parabolic cylinder function form D„(|g|) with 
q = (xi — X2)/x osc [18j|. The argument establishing the 
connection with the hard sphere gas in this case is the 
same as given previously for spinless bosons [l7|- The 
only difference between the sTG and HSB wave functions 
apart from normalization is that the sTG wave function 
allows penetration into the region interior \x\ < dh.s.- 
This holds for — oo < A < 0, and the penetration is small 
when |A| 3> 1 so that x no d e = ih.s. ~ ^id- Furthermore, 
the sTG and hard sphere solutions have the same energy 
since they satisfy the same Schrodinger equation in the 
exterior region, and they have the same en ergy eig envalue 
provided that the energy index v of D v [lj, LL8| is cho- 
sen to be the same as that for the ground state solution 
D v in the entire region \x\ > including the attractive 
delta function at the origin. The sTG wave function D u 
is highly excited relative to the collapsed McGuire state 
due to the effects of the strong attraction at x = 0, but 
when restricted to the region |x| > Xnodc, it is identical 
with the ground state of the hard sphere Bose gas apart 
from normalization, and the energies of the highly excited 
sTG state and the hard sphere ground state are exactly 
equal HI. 

A generalization of this theorem should hold for all 
N > 2. Consider a generalized hard sphere system with 
a hard core interaction of diameter a; no d e for opposite spin 
pairs, but no interaction for (ft) and (Li) pairs. For op- 
posite spin pairs, the ground state wave function of this 
system will vanish when \xj — X(\ < x no( jc but be nodeless 
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for \xj — X(\ > Xnodo, whereas for (ft) and (4-4-) pairs it 
will be of the trapped ideal gas form. It therefore seems 
likely that the unknown exact sTG energy eigenstate, to 
which (|10p is an approximation, will differ from such a 
generalized hard sphere ground state only in the interior 
region \xj — xt\ < x no de for opposite spin pairs, and in 
overall normalization. Neither the sTG state nor the gen- 
eralized hard sphere state are known exactly. Therefore I 
state the belief of the identity of the sTG and generalized 
hard sphere wave functions and energies for N > 2 as a 
conjecture. 

There also exists a quite different sTG state when 
N+ = iVj. = N/2. In this case it was shown recently 
12], for a system on a ring with no circumferential trap- 
ping potential, that for — oo < A <C — 1 there exists a 
gas-like energy eigenstate in which all the fermions are 
tightly bound into (ji) dimers which behave as strongly 
attracting bosons, forming an exact analog of the un- 
trapped solution for spinless bosons [24j. It is a gas-like 
sTG state lying much lower than the fermionic sTG state 
of Eq. ([TO"]) due to the negative binding energy of the 
dimers, but much higher than a McGuire-ideal Fermi hy- 
brid state similar to that of Eq. © . One expects that a 
similar sTG state exists for the present case of harmonic 
trapping. In fact, for N = 2 it has already been given as 



the JVf = Ni — 1 case of the previously-discussed exact 
N = 2 solution, and for larger even N there will be a 
Bose-like sTG state similar to that of (24j. However, it 
could not be created by a sudden switch from a strongly 
repulsive to strongly attractive interaction as in (16J. In 
fact, the sudden approximation probability of finding the 
system in such a state after a switch A^> 1 — >• — A^; — 1 
is exactly zero in the limit |A| — ¥ oo, since the state (|10[) 
is an exact energy eigenstate for both A = +oo and 
A = — oo. Then by continuity this probability will be 
<C 1 after a sudden switch A 3> 1 — > —A <C — 1. 
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